Spontaneous creation of discrete breathers in Josephson arrays 
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We report on the experimental generation of discrete breather states (intrinsic localized modes) 
in frustrated Josephson arrays. Our experiments indicate the formation of discrete breathers during 
the transition from the static to the dynamic (whirling) system state, induced by a uniform external 
current. Moreover, spatially extended resonant states, driven by a uniform current, are observed to 
evolve into localized breather states. Experiments were performed on single Josephson plaquettes 
as well as open-ended Josephson ladders with 10 and 20 cells. We interpret the breather formation 
as the result of the penetration of vortices into the system. 
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Discrete breathers (DBs) are spatially localized, time- 
periodic excitations of spatially uniform nonlinear lat- 
tices. After initial theoretical studiesii^, they were ob- 
served in experiments in a variety of systemsii*^. In 
arrays of superconducting Josephson junctions, so-called 
Josephson ladders (JL), DBs can be excited in a con- 
trolled vf&y^i^^. Experiments carried out in JLs ranging 
from one cell up to several tens of cells revealed a large 
variety of possible DB states^, studied their region of ex- 
istence and the resonant interaction of the nonlinear DB 
excitation with small-amplitude linear wavea^iifl. 

DB states in JL are usually excited from the uniform 
static state of the lattice by using additional currents 
(i.e., local forces) which are removed after the DB cre- 
ation. This well-controlled procedure helps to study DB 
properties systematically. In this paper we report that 
DB can be also excited without applying any local forces 
to the lattice. We present several scenarios where DB 
states are observed to form spontaneously. The DB for- 
mation may occur either from the initially static system 
state, or as a "breakdown" of spatially extended whirling 
states that are in resonance with the cavity eigenfre- 
quency. 

Within the RCSJ (resistively and capacitively shunted 
junction) model^^ a Josephson junction (JJ) has a simple 
mechanical analog, which is a pendulum. In this analog, 
the difference of the superconducting phases across the 
JJ corresponds to the tilting angle of a pendulum. The 
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FIG. 1: Sketch of the single cell system (a) and of the Joseph- 
son ladder (b). Josephson junctions are indicated by crosses, 
superconducting leads by straight lines. The homogeneous 
bias current, /, is injected through resistors Rb- 



current sent across the JJ is represented by a torque 7 
applied to the pendulum. Both JJ and pendulum are 
subject to a velocity dependent friction force aip with 
damping constant a. Finally, the JJ or pendulum phase 
moves in a cosine potential due to the Josephson energy 
or gravitation, respectively, which leads to a force pro- 
portional to siTup. The dynamical equation governing 
the phase motion can in either case be written in nor- 
malized units as 



ip + atp + sm (/3 = 7. 



(1) 

The (approximate) time-averaged solutions {^p) and {lp) 
of this equation may have several simple forms. For 7 > 
1, Eq. Q has only a whirling solution with {ip) — 7/ a <^ 
if) — (7/0^)^- For the intermediate range 4a/7r < 7 < 1, 
a static solution with [ip) = arcsin7 coexists with the 
whirling solution. For 7 < Aa/i:^ only the static solu- 
tion survives. Small-amplitude oscillations of may be 
considered as perturbations to the time-averaged solu- 
tions. Due to the Josephson law, the time derivative of 
the phase is proportional to the voltage drop across the 
JJ. For the whirling solution, the normalized voltage v is 
therefore proportional to the current, v = (ip) = (l/a)7, 
and the inverse of the damping a is simply the (nor- 
malized) ohmic resistance r of the J J, r — 1/a. The real 
current / is normalized to the maximum current allowing 
for the static state Ic, hence 7 = I/Ic, for convenience. 

A JL is a single row of a two-dimensional rectangu- 
lar Josephson junction array (see Fig. ^). The JJs are 
connected via superconducting leads. JJs on the spars of 
the ladder will be referred to as the "vertical" junctions, 
those along the rungs as "horizontal" junctions. The ver- 
tical JJs are all identical and have an area of A^j , likewise 
all horizontal JJs are identical but have an area Ah- A 
single cell of the ladder consists of two vertical and two 
horizontal JJs. Many of the properties from the extended 
Josephson ladder are preserved in a single cell with only 
one horizontal J J (Fig.QJi). This system is the simplest 
possible Josephson array permitting DB-like excitations. 
The anisotropy 77 is defined as the ratio of the horizontal 
and vertical J J areas, 77 = A^/Ay. The J J phases inside 
a cell are coupled due to the magnetic flux quantization 
and via a circulating inductive current. Every cell has a 
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self-inductance L. A simple model"'^- describing the phase 
dynamics of the JL neglects mutual inductances between 
cells and is similar to the discrete sine-Gordon model. 

The vertical J Js form an array of rotors (pendulums) . 
These JJs can be driven by a uniform bias current. Tech- 
nically, this is achieved by using a single current source 
which feeds a comb of parallel resistors Rb (see Fig^). 
If the resistance Rb is large enough, the currents through 
all vertical branches are identical. If the normalized cur- 
rent 7 per vertical JJ is large compared to unity, all 
vertical JJs rotate. This is the so-called homogeneous 
whirling state (HWS). In contrast, if the current is small, 
all JJs remain in the static state (SS). In the region 
4a/7r < 7 < 1, the ladder may have some vertical JJs 
in the whirling state, while the remaining JJs are in the 
static state. In the standard terminology, this localized 
state is called roto-breather state (RB). Because of to 
the JL geometry, localized voltage states always have to 
conform to Kirchhoff's voltage law. Hence a single ver- 
tical rotating JJ has to be accompanied by one or more 
horizontal JJs in the whirling state. A large variety of 
different RB states is possible, as it has been observed in 
experimenliSii^. 

The SS of the ladder is influenced by an externally 
applied magnetic flux $ext = Bcxt^ (due to a magnetic 
field -Bext threading the cell area A) . The magnetic field 
leads to a screening (Meissner) current flowing around 
the JL. The field may penetrate into the ladder, lead- 
ing to vortices in the screening current distribution. A 
measurement of the maximum critical current 7c of the 
S state reveals information on the magnetic field pene- 
tration into the ladder. In Ref. I 1. the existence region 
of the S state of a Josephson ladder was investigated and 
compared to theoretical predictionsi^. Deviations were 
found for close to half-integer frustration and attributed 
to the presence of meta-stable states. In this paper, we 
present similar measurements of 7c vs. magnetic field 
and will argue that they are related to the creation of 
DBs. 

When the ladder is in the HWS, all vertical J J phases 
rotate at an identical average frequency fJ, while the hor- 
izontal JJs do not rotate. This is accompanied by a uni- 
form voltage drop across the vertical JJs. Superimposed 
phase oscillations are strongly enhanced if the rotation 
frequency $7 matches one eigenfrequency of the cavity 
formed by the finite ensemble of cells. The enhanced 
oscillations lead also to an enhanced damping. A suit- 
able experiment is to measure the differential resistance 
(which is inversely proportional to the damping) of the 
JL. It is defined as the ratio of voltage change induced 
by current variation. The matching frequencies (volt- 
age positions) are characterized by steps in the current 
vs. voltage characteristics. Such resonances of the HWS 
with cavity modes have been studied by Caputo et ali^. 
A similar behavior was observed when frequency compo- 
nents of a localized RB state coincide with one of the cav- 
ity eigenfrequencies®'^^. We will show experiments where 
the HWS resonances are found to enforce the creation of 



symmetry-broken RB states. 

We performed the first group of experiments on a sin- 
gle cell, having the central hole area A — i x 4 ^m^, 
with three small underdamped Nb/Al-A10a;/Nb JJsii, 
Fig. Hi. The vertical JJs have an area A^ = 6 x 6 /im^, 
while the horizontal J J has an area Ah — 3 x 4 /im^. 
The anisotropy of the system is then rj = Ah/ Ay — 0.33. 
At 4.2 K the parameter of self-inductance evaluated by 
measurement of a SQUID of the same geometry, is (3l = 
2ttLIcv/^o = 0.98. Measurements were done in a tem- 
perature range between 4.2 K and 6.3 K. At each temper- 
ature position the damping a = 7r7r/4 and the critical 
current density jc are evaluated from the current-voltage 
characteristic of HWS. At T = 4.2 K, a ~ 0.03 and 
jc — 126 A/cm^. The bias current / is introduced via 
two parallel external resistors, Rb = 1.5 kfi, see Fig. ^ 
so that the bias currents on the two vertical branches are 
equal. For convenience, the experimental data are shown 
in the normalized units of the bias current 7 = I / lev and 
frustration / = ^cxt/^o, where $0 is the fiux quantum. 

We ramped up the bias current 7 and registered a tran- 
sition of the system to a resistive state by detecting a 
non-zero voltage at either junction Vi (this critical value 
of 7 we call jI) or junction V2 (we call it 7c). The 
measured dependence of 7c and 7c on the external mag- 
netic field is, as expectedii, periodic, see Fig. |2| Around 
/ = ±0.5, the behavior of the two vertical JJs is found 
to differ and spontaneous formation of broken-symmetry 
statesiS is observed. In Fig. |21 the experimental data 
at r ~ 5.5 K are shown. In the range of frustration 
-0.7 < / < -0.5 and 0.3 </ < 0.5, where large clock- 
wise mesh currents are induced, the junction Vi remains 
in the SS while the junction V2 switches to the resistive 
state. The junction Vi switches to the resistive state at a 
higher current about 7 ~ 0.8 that limits the existence of 
this broken-symmetry breather stateiSii^. As expected, 
the maximum value of the bias current that allows for 
the existence of the breather state does not depend on 
the frustratioriiS. In the range —0.5 ^ / ^ —0.3 and 
0-5 ^ / ^ 0.7, large counterclockwise mesh currents are 
induced. The junction Vi switches first and the second 
possible breather state is spontaneously induced. The 
range of frustration, where the magnetic field induced 
breather states are observed, broadens with temperature. 

At temperatures higher than 4.2 K, the HWS current- 
voltage characteristic in presence of external frustration 
showed the well-known resonance due to the interac- 
tion with the electromagnetic oscillatory modes of the 
celliS*i2iS£. At T ~ 5.5 K, for some values of frustration 
about 0.5, the interaction of the HWS with the oscilla- 
tory modes of the cell led the system to switch from the 
top of the resonance towards a RB state, see Fig. 

The second group of experiments has been performed 
with open-ended JLs. Measurements of the critical cur- 
rent needed to force the vertical JJ i to the finite 
voltage state are presented in Fig. ^ The sample con- 
sisted of 20 cells, with parameters rj = 0.5, Pl ~ 1-2 and 
a « 0.03. The sample layout is similar to the one pre- 
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FIG. 2: Dependence of the depinning currents jl and 7^ on 
the external frustration, shown respectively by full and open 
circles. The hatched patterns show the regions where junc- 
tions Vi and V2 are in the SS. 



FIG. 4: Dependence of the left (open squares, 7^), central 
(filled diamonds, 7^^), and right (open circles, 7^^) vertical 
junction depinning currents on the frustration /. 




voltage on V^, (mV) 



FIG. 3: Dynamic switching from the resonance of the HWS 
(black dots, with a resonant step at about 0.2 mV) to the 
broken symmetry state (gray points); / = —0.41 



sented in Ref. 14. The measurement was performed at a 
temperature of 4.2 K. 

The pattern of the depinning current 7* is periodic in 
frustration, with the period corresponding to a magnetic 
flux of $0 per cell. Around integer values of /, the critical 
currents for left, center, and right positions in the ladder 
coincide. After the depinning, for 7 > 7J, the HWS is 
formed. For / 0.5, the critical currents 7* differ for 
voltages measured on different positions in the ladder. 
At —0.35 < / < —0.25, the right and central vertical JJs 
of the ladder start to rotate at a lower value of 7 than 
the left part of the ladder. For -0.65 < / < -0.35, the 



central part depins at the lowest value of the current, the 
right part at an intermediate current, and the left part at 
the highest value of the current. For —0.75 < / < —0.65 
the center and left parts depin for a low current and the 
right part depins at a higher current. 

By definition, a breather state is formed when one part 
of the vertical JJs in the ladder rotates while the rest of 
it remains in the SS. Hence in the described frustration 
regions around / ^ 0.5, breathers are found. They form 
solely due to an external magnetic field. Because of the 
limited number of voltage probes in our long ladder, we 
were not able to exactly identify the precise number of 
rotating JJs in the breather states which were formed, 
but could clearly distinguish them from the HWS. 

Our observations can be qualitatively explained in the 
following way. When the HWS has formed for large 7 and 
finite frustration /, magnetic fiux penetrates into the sys- 
tem via the boundary cells. When 7 is reduced to zero, 
the ladder settles in the SS, which may be a meta-stable 
state or the system ground state. Above some critical 
value of the frustration /, this state consists of a nonzero 
amount of vorticesi^. For small bias currents 7, the vor- 
tices are pinned due to the lattice discreteness. However, 
when the vortices start to move, they may ultimately 
bump into the ladder boundary. This will eventually 
force the boundary JJ to the whirling state, with the rest 
of the ladder remaining in the static state. The low value 
of 77 is crucial for this behavior. Preliminary numerical 
simulations have shown this kind of behavior for typical 
experimental parametersSi. For longer systems, it is also 
expectable and observed that moving vortices may be ex- 
pelled via one of the horizontal JJs before reaching the 
boundary. When a vortex passes "through" a horizontal 
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voltage, V(mV) 



FIG. 5: Current-voltage curve of a 10-cell JL indicating the 
creation of a breather state from a whirling mode resonance 
for a frustration / = 0.66. The current was swept from 
1 = 8 fiA to 7 = 25 and back. Both the central (circles) 
and the right boundary (diamonds) voltages are plotted. Ar- 
rows indicate directions of the voltage switching and the path 
which is taken during the current sweep. The square indicates 
the initially prepared resonant HWS (/ = 8.4 pA). The circle 
indicates the position at which the symmetry-broken state is 
formed (7 = 21.4 fJ,A). The inset shows the dependence of the 
center junction depinning current on frustration, the arrow 
marks the frustration at which the current-voltage curve was 
taken. 



junction, this phase may be left rotating, again leading 
to the formation of a RB. A series of DBs may be left be- 
hind a moving vortex. This has been observed also in the 
simulations presented in Ref. [2^ We conclude that the 
depinning transitions for JLs with low anisotropy rj can 
lead to the formation of localized whirling (RB) states. 

We also observed the creation of a RB state from a 
resonant HWS in JLs. These experiments were per- 
formed with a 10 cell open-ended JL of anisotropy rj = 
0.49. Other parameters during the measurement were 
Pl ~ 0.35 and a ~ 0.03. In the presence of the ex- 
ternal magnetic field, we observe resonant steps on the 
current-voltage characteristics of the ladder, as shown in 
Fig. |31 They are occur when the vertical J J rotation fre- 



quency n matches one of the cavity eigen-frequencies, as 
discussed in Ref. When the homogeneous current is 
increased at a resonance, the amplitude of the cavity os- 
cillations increases. If the driving current becomes too 
large, the resonant state cannot be maintained anymore 
and the non-resonant HWS (now with a larger rotation 
frequency) is retained. This hysteretic process is usually 
repeatable for an arbitrary number of times. However, 
in some rare but reproducible cases, we observed a tran- 
sition from such a resonant, spatially homogeneous state 
towards a localized RB state. The localized state RB 
can be identified as before by a zero voltage drop across 
some vertical JJs in parallel with a finite voltage drop 
across another group of vertical JJs. In Fig. |S1 the volt- 
age switching, and hence the RB formation happens in 
the region emphasized by the circle, around / = 21.4 fiA. 
During further increase of the current the RB voltage 
stays constant, apart from some tiny jumps. We thus 
observe here that a resonant RB state is formed from 
the resonant HWS. Ultimately, the RB voltage switches 
to a high value (at / « 25 /lA). Here the breather state is 
still present, but is non-resonant. The tiny voltage jumps 
along the RB resonance may indicate either transitions 
between different excited modes or a change in the RB 
size, which could not be detected in this experiment due 
to the limited number of voltage probes. The described 
transition occurs for a very narrow range of magnetic 
field. 



In summary, we reported here on the spontaneous 
creation of discrete breather states in single Josephson 
plaquettes and Josephson ladders. The excitation of 
breather states is observed during the transition from the 
static to the whirling state of the array, resulting from 
a ramp of an externally applied homogeneous bias cur- 
rent. In another regime, the breather states are created 
from the resonant uniformly whirling state of the array. 
The spontaneous creation appears in the presence of an 
external magnetic field close to half-integer frustration, 
which leads to large screening currents. The observed 
behavior may be related to the appearance of symmetry- 
broken states in two-dimensional Josephson arraysi^^. An 
analysis of the described phenomena could provide more 
insight in the physical relevance of localized voltage states 
in Josephson coupled systems. 
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